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Abstract The technique of projecting the four-dimensional two-body Bethe-Salpeter 
equation onto the three-dimensional Light-Front hypersurface, combined with the quasi- 
potential approach, is briefly illustrated, by placing a particular emphasis on the re- 
lation between the projection method and the effective dynamics of the valence com- 
ponent of the Light-Front wave function. Some details on how to construct the Fock 
expansion of both i) the Light-Front effective interaction and ii) the electromagnetic 
current operator, satisfying the proper Ward-Takahashi identity, will be presented, ad- 
dressing the relevance of the Fock content in the operators living onto the Light-Front 
hypersurface. Finally, the generalization of the formalism to the three-particle case will 
be outlined. 

Keywords Light-Front Field Theory • Bethe-Salpeter equation • Few-body systems • 
Electromagnetic structure 



1 Introduction 

In his 1949 seminal paper [T], P. A.M. Dirac proposed three peculiar representations of 
the Poincare group, in strict relation with the choice of possible space-time hypersur- 
faces without a time-like direction, (see also, e.g., [2]): the Instant form, the Light-Front 
(LF) form and the Point form (cf [3] for a recent review). Each hypersurface leads to 
a specific, but equivalent, description of the dynamics of a relativistic interacting sys- 
tem. As a matter of fact, the choice of an initial hypersurface, in which the points 
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are separated by space-like distances and therefore no causal connections are allowed, 
suggests the set of dynamically independent variables, suited for implementing a de- 
scription of the initial state of an interacting system. In conclusion, Dirac indicated 
the potentially fruitful role of the Hamiltonian approach for relativistically describing 
interacting systems, within a field theoretical framework, as well (see, e.g., [2]). 

In the Instant form, corresponding to the choice of a constant-time hypersurface 
in the Minkowski space (xq = 0), translations and rotations, being the generators 
in the stability set of the chosen hypersurface, commute with the Hamiltonian. The 
eigenstates of the Hamiltonian, constructed in terms of the Fock-space basis (in this 
basis the free Hamiltonian is diagonal), have both the eigenvalues of the total three- 
momentum and the total angular moment as good quantum numbers. Translations and 
rotations do not mix different Fock-components of the wave function and therefore a 
truncation in the Fock-space, that is required for practical applications, is stable under 
those transformations. Differently, the Instant-form boosts have dynamical nature and 
therefore they mix the Fock components. The shortcoming of a truncated Instant-form 
Fock basis is made manifest by the lack of Lorentz invariance, e.g. when we evaluate 
the expectation value of observables involving initial and final rest-frames, having a 
relative non zero velocity. 

In the LF form, where the hypersurface x + = t + z = is chosen for quantizing 
the theory, one has a kinematical subgroup of the Poincare group, built by seven 
generators (see e.g. [21S1E1E] ) - Such operators, keeping invariant the LF hypersurface, 
do not contain the interaction and they correspond to i) three LF translations, ii) 
a longitudinal LF boost along the z-direction, iii) a rotation around z, and iv) two 
transverse LF boosts, suitable linear combinations of transverse Instant-form boosts 
and rotations. The remaining three generators have a dynamical nature: i) two LF 
transverse rotations, like the LF boosts but with different sign combinations, and ii) 
the Hamiltonian P~~ = P° + P 3 , i.e., the generator of LF-time x + translations. From 
this classification, immediately an important feature stems out: the stability under 
LF boosts of any truncated Fock expansion of the physical states (eigenstates of the 
squared mass Fock-space operator, see, e.g., 0), since the LF boosts are diagonal in the 
Fock space [7]. This property is fundamental for a consistent treatment of the boosts 
when a truncated Fock basis is adopted. For example, if one calculates electromagnetic 
form factors for momentum transfers along the z— direction ( q + > and = [5]), 
the initial and final rest-frames are related simply by a kinematical transformation, 
and therefore one could use a truncated description of the initial and final states. As a 
last remark, it is worth noting that the change to LF variables allows one to linearize 
the dependence upon the dynamical variable (conjugated to the LF time), with a 
simplification in the analysis of the poles in the propagators (i.e., fcg — (m 2 + |k| 2 ) + ie — > 
k + [AT - (m 2 + \kj_\ 2 )/k + ]+ie). 

The physical outcome of a quantum theory does not depend on the particular 
space-time hypersurface adopted for quantizing, and therefore the theory has to have 
an equivalent four-dimensional (4D) formulation, where explicit covariance is built 
in. In this respect, the Bethe-Salpeter equation (BSE) is a 4D field-theoretical tool 
[5] to explore the non perturbative physics of bound and scattering states of few- 
particle systems. Within a LF framework, it can be easily shown that the projection 
onto the LF hyperplane of the BS amplitude is proportional to the valence component 
of the physical state (see, e.g., [5]). This suggests to take the BSE as the starting 
point for an investigation of the internal dynamics, alternative to the one where a 
coupled-equation system, generated by the Hamiltonian acting on the full Fock space, 
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is considered. In order to pave the way from the BSE to the eigenequation for the 
3D valence component and back, another important suggestion is offered by the idea 
of "iterated resolvents" of Refs. |1UU11| . where it is proposed that the complexity of 
the Fock-space LF Hamiltonian builds an effective dynamics for the 3D valence state. 
Indeed, it should be pointed out that the mass operator for the 3D valence component 
can be exactly obtained from the BSE by a quasi-potential (QP) expansion [12], as 
investigated in Refs. [13U14] and [TS] (where the associated set of coupled equations for 
LF Green's functions was derived). 

The QP approach makes feasible to single out the "trivial" global propagation 
of the interacting system by means of a well-defined auxiliary Green's function (see 
Eqs. ([9]) and (JT3J) , below). This fundamental goal, in turn, allows one to reconstruct 
the BS amplitude from the 3D valence component, with a one-to-one correspondence, 
and moreover, to obtain the expansion of the 3D effective interaction, entering the 
eigenequation for the valence component. Each contribution to the expansion immedi- 
ately acquaints a transparent physical interpretation (given the strong analogy between 
the LF evolution and the non relativistic one) that produces a straightforward guidance 
in evaluating the relative importance of the associated diagrams. The same analysis 
seems more involved in the studies of the Hamiltonian in a truncated Fock-space, i.e. in 
the analysis of the generated system of coupled equations. Noteworthy, the convergence 
of the expansion, that has a direct impact on the lost of covariance with respect to 
the subset of the dynamical transformations, has been investigated in a simple bosonic 
Yukawa model in Ref. [13] . where the Fock content of each contribution has been recog- 
nized as the ordering "parameter". In closing, the QP technique appears so appealing 
within the LF framework, that it could be very interesting to implement comparisons 
between actual calculations performed in other relativistic QP approaches, like the Co- 
variant Spectator Theory (see e.g. [22]) applied in many relevant few-nucleon problems. 

The LF projection of the BSE for few-particle systems, and the consequent trunca- 
tion of the physical state in the Fock space, can be useful if there is a dominant valence 
state, or if the normalization can be saturated at large extent by including only few 
Fock components. In Nuclear Physics, the nucleonic component is largely dominant, 
as in the deuteron case, while in applications to hadrons, Fock components beyond 
the valence one have been recognized to be relevant, in particular in the description of 
physical quantities pertaining to inelastic channels, like structure functions and gen- 
eralized parton distributions (from deeply virtual Compton scattering [ID]). Recently, 
a signature of the the relevance of the components beyond the valence one has been 
singled out in the evaluation of the nucleon electromagnetic form factors by using con- 
stituent quark degrees of freedom in the = frame [5]. As a matter of fact, a zero in 
the proton electric form factor has been associated with a cancellation between valence 
and nonvalence contributions to the electromagnetic current (see e.g. [17] and [18] for 
the experimental overview). 

Aim of the present review is to yield some insights on i) the technique for projecting 
the 4D BSE onto the three-dimensional (3D) LF hypersurface using the quasi-potential 
approach and ii) the construction of the effective dynamics of the valence component of 
the LF wave function for both two- [131114] and three-particle systems [19] ■ Moreover, 
some details on how to obtain an effective electromagnetic current operator i) acting 
on the valence component of composite two-boson and two-fermion systems, and ii) 
fulfilling the Ward-Takahashi identity (WTI) [2D] and [31], are given. 

The review is organized as follows. In Sec. 2 the QP reduction is presented as a 
tool to eliminate the relative LF-time. In Sec. 3, the relation one-to-one between the 
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BS amplitude and the LF valence wave-function is illustrated. In Sec. 4, a hierarchy 
of coupled equations for LF Green's functions is discussed in order to show the Fock 
content of the projection of the BS equation onto the LF hyperplane. In Sec. 5 the 
projection technique for three-body BS equations is briefly discussed, along with the 
role of induced three-body forces. In Sect. 6, the electromagnetic current operator, 
acting on the valence component and fulfilling the LF Ward-Takahashi identity is 
introduced, with some remarks on both structure and covariance of the truncated (in 
the Fock space) current. In Sec. 7, a summary and some perspectives are presented. 



2 The Quasi-Potential reduction: a tool for eliminating the relative 
LF-time 

The 4D BSE (see, e.g. [9]) is the starting point of many studies of few-body systems 
which aim to account for relativity. The BSE requires a relativistic field-theoretical 
approach, based on an interacting Lagrangian able to model the system under investi- 
gation. In the particular case of a two-body scattering, the 4D scattering equation for 
the transition matrix, T(K), with total four-momentum K is written as follows 

T{K) = V(K) + V(K)g (K)T(K) , (1) 

where the interaction V(K) contains, in principle, all the possible two-body irreducible 
diagrams. The two-body disconnected Green's function, Qq(K), should include self- 
energy terms, but they are neglected in the approach developed so far. Therefore, in 
the case of two bosons, 9o(K) becomes 

i 2 1 1 
G (K) = ^ , (2) 

^ n k\ — mf + is k\ — m?, + is 

where k^ is the four-momentum operator and the factor 2-k is introduced for conve- 
nience. For two fermions one has Qq(K) — ► Gq(K) = f^j + mij (jt 2 + m 2^ Gq(K). 

The two-particle bound-state with total 4-momentum Kg, K b = Mg, corresponds 
to a T-matrix pole. The residue is associated with the vertex function, namely the non- 
trivial of part the BS amplitude, tyg , solution of the following homogeneous equation 

Wb) = Go(Kb)V(K b ) \* b ) . (3) 

The normalization condition [9] has to be satisfied in order to fully determine \^b) ■ 
For scattering states with total 4-momentum K, the BS amplitude is a solution of the 
following inhomogeneous equation 

|^ + ) = |#o) +Gq(,K)V{K)\1> + ) , (4) 



while the T-matrix, solution of Eq.([T]), corresponds to the connected four-point function 
which brings information on both the scattering and bound states. In Eqs. and @ 
the four-momentum conserving 5-function is factorized out. 
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Since in the LF projection method a central role is played by the on-minus-shell 
propagation, as it will be emphasized below, let us write the relevant matrix elements 
of the free two-boson Green's function, viz 



(k[-\G (K)\k^) = 

i 2 8 (k[- - k- 



271 k+(K+ -k+)[k-- k- on + ^)[ K --k-- k- on + 



k+ / \ 1 2on K+-k+ 

,„3 u + — 1,0 , ;„3 



(5) 



i± ) 



where the LF four-momenta are k^ = (k~ := fe" — kf , k^ := k^ + kf , k 
k~ on = (k;j_ +m 2 )/k^ (i = 1, 2) is the on-minus-shell momentum operator, with eigen- 
functions given by the LF plane waves, (x~n i± |fc+k i± ) = Af e~ l ^ k t x i - k ii- x i-L). 
The completeness relation and the normalization are 

and (k' + k' ± \k + k ± ) = 2(2ir) 3 S(k' + - k+)5{k' ± - kj_), respectively. 

The free two-fermion propag ator, Gq (K), can be decomposed in an on-minus-shell 
term, Gq(K), and a part that contains the so-called instantaneous (in the LF-time) 
contribution, since the Dirac propagator can be separated in two terms as follows 

]/• _1_ Tii ]/• „L_ in "i/ - ^ - 



k 2 - m 2 + is k+(k- - ko„ + ^) 2fc+ 

where the first term yields the on-minus-shell propagation, while the second one the 
LF-time instantaneous term of the Dirac propagator. For the fermion case, we will be 
interested in the analogous of Eq. ([5]), written in terms of Gq(K), namely 

(k[- 1 com |fcr) = (j lon + mi ) (£ 2on + m 2 ) (k'r \ g (k) i^-) , w 

For the sake of a unified formal treatment of two-boson and two-fermion systems, in 
■what follows we put Gq(K) = Gq(K). 

One could easily extend the present analysis to systems composed by particle- 
antiparticle or by other mixtures, like a fermion and a boson (see e.g. [23]). 

The first step for projecting the BSE onto the LF surface is the introduction of the 
free-resolvent, i.e., the Fourier transform in K~ of the global a; + -time free propagator 
of the two-particle system. This amounts to integrate the matrix elements, Eq. (|5j) (or 
Eq. (©), of the 4D Gq(K) over fc^~ and so that the relative LF time between the 
particles is eliminated, and one remains with a dependence upon K~ , i.e. 

\G (K)\ := J dk'fdk^ (k'f\G (K) |*i ) = g (K) (9) 

where go(K), called the free-global LF propagator, is a 3D operator depending upon 
the LF momenta (fc^~,kjj_) only, and it is explicitly given by 

9o(K) = — — ^ ^ , (10) 

k+(K+ - k+) (K- - k- on - k- on + l e) 
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where P = i9(K + — k^)6(k^) for two-bosons and 

P = i6{K+ -k+)e(k+)(2m 1 )(2m 2 )A + (k lon )A+(k 2on ) 

for two-fermions. The projector A+(k on ) = (jt on + m\ /2m is the positive energy 

spinor projector. A positive value for K + is used without any loss of generality. 

In Eq. |[9j), the vertical bars "| "on the right side and on the left one indicate that 
the minus components in \k~) and {k'~ | have to be integrated out [13II14] . respectively. 
Notice that, for two fermions, the inverse of go(K) exists only in the valence sector, 
since the projectors, single out only positive energy states. 

Within the QP approach [12], where an auxiliary interaction, W(K), is introduced, 
the full T-matrix is solution of the following coupled equations 

T(K) = W(K) + W(K)G (K)T(K), (11) 
W(K) = V(K) + V(K)A (K)W(K) , (12) 

where A Q (K) = G (K) - G (K) for bosons and A (K) = Gq {K) - G (K) for 
fermions. The auxiliary Green's function Gq(K) is a 4D operator, depending upon 
the four- momenta of the two constituents, and it represents the key quantity of the LF 
projection. It is the 4D image of the 3D dimensional go(K), that, we strongly stress, 
does not contain the relative-time propagation of the system. Therefore Aq(K) just 
takes into account such a propagation in the 4D space, namely it will be an essential 
ingredient in the description of the internal dynamics of the systems. The 4D operator 
Gq(K) is defined by 

G (K) = n (K)go(K)n (K) , (13) 

where 

II (K) = G (K)\ 9o\K) , n (K) = g^(K) \G (K). (14) 

These operators, the free reverse LF -•projection operators, connect three and four di- 
mensional quantities. In the next section the corresponding interacting operators will 
be introduced. It is worth noting that the choice of Gq(K) and the corresponding 
integral equation for W(K), Eq. ()12|) . allows only for LF two-body irreducible terms. 
The solution by iteration of Eq. (|12|l is given by 

oo 

W(K) = ^ W n (K), (15) 

n=l 

with W(K) n = V(K) (A (K) ViK))^ 1 . The diagrammatic analysis of the series ([To]) 
shows that for each term one has a specific Fock content, associated to the propagation 
of the virtual intermediate states, as discussed to some extent in Sect. [4] Of course, 
a truncation of the sum in Eq. (|15[) puts a bound on the number of the Fock compo- 
nents involved in the actual calculation. Moreover, the same holds for the 3D effective 
interaction w(K) := IIq(K)W (K)IIq(K) that determines the valence component, as 
shown in the following Section. In particular, the 3D effective interaction, since it is 
non diagonal in the Fock space, makes possible the coupling of the valence sector to 
the higher Fock components of the wave function (formal details are given in Sec. 4 
and [IS]). 
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It is worth noting that, in model studies (sec 13 , 24 , 25 , 26 , 27, 28 ), the expansion 
(]15[) is rapidly converging, since the probability of higher Fock states is quickly de- 
creasing [2U1I5D] . 

In fermionic models, with point-like couplings, the LF BSE was investigated by 
retaining the lowest order kernel and subtle problems, related to the divergences pro- 
duced by the dependence upon the transverse momentum, were found |31H32ll33] . Part 
of the difficulties can be ascribed to the absence of the instantaneous terms as analyzed 
in Refs. (141133] . 

Following Refs. [I3ll 14U20U21 j . one can construct a 3D LF T-matrix, t(K), from the 
4D one. In particular, one has 

t(K) = n (K)T(K)n (K) = w(K) + w(K)g (K)t(K) = w(K) + w{K)g{K)w(K) 

(16) 

where g(K) is the interacting LF Green's function, fulfilling the integral equation 

g{K) = g (K) + g (K)w(K)g(K), (17) 

Notice that it also holds g(k) — go(k) + go(k)t(K)go(k). The 3D operator, g(K), is the 
Fourier transform in K~ of the global LF-time propagator, viz 

g(K) = \G (K)\ + \G (K)T(K)G (K)\, (18) 

and evolves the system from an initial state, defined on a given LF hypersurface, to 
another one, after a LF-time interval xt — xf > 0. By iterating once the integral 
equation (|lip . and using Eqs.fTJ]) and (|16p . one has 

T(K) = W(K) + W(K) G (K) + G (K)T(K)G (K) W{K) = 

= W{K) + W(K)7l (K)g(K)n (K)W(K) , (19) 

This relation allows one to map the 3D dynamics into the 4D space. 

It turns out that the on- mass-shell matrix elements of T(K), which define the 
two-constituent scattering amplitude, are identical to the ones obtained from the on- 
minus-energy-shell matrix elements of t(K) (see discussion in [1311 14] ). Unless otherwise 
indicated, the operators go(K) and w(K) have to be evaluated with a "+ie" prescrip- 
tion. 



3 The BS amplitude and the LF valence component 

The relation between the BS amplitude \&b) an d the valence component of the LF 
wave function, for a bound state with total momentum Kg is given by Q3] (let 

us recall that we formally put Gq{K) = Gq(K)) 

Wb) = Go(K b )W(K b )G (K b )\ go\K B ) \4> B ) , (20) 

where the valence component is the solution of the following eigenequation 

\$b) =9o(K B )w(K b )\$ b ) , (21) 

It should be pointed out that, in the case of fermions, the instantaneous terms from 
the Dirac propagators appear in i) Qq = Gq , ii) W and iii) the effective interaction w. 
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The identity G (K B )\ (gQ 1 (K B ) - w{K B )) \<P B ) = can be added to Eq. (20J in 
order to get 

\$ B ) = [1 + A (K B )W(K B )] G (K B )\ go\K B ) \$ B ) , (22) 

This expression holds not only for bound states but also for scattering states (see e.g. 
HB), viz 

= n(K) |<Z>+) (23) 
where 17 (K), the interacting reverse LF-projection operator, is given by 

n(K) := [1 + A (K)W(K)] G (K)\ g^{K) = G R (K)\ g~\K) . (24) 

with 

G R (K) := G (K) + G(K)V(K)G (K) = G(K)G^ (K)G (K) , (25) 

Notice that in G B (K) the on-minus-shell-Green's function Gq(K) appears on the right- 
most position, and this leads to apply the "|" operation on the right. 

The operator IJ(K) acts on the valence component of the LF wave function in 
Eqs. (|22|) and (|23[) and it allows one to fully reconstruct the 4D BS amplitude for 
both bound and scattering states, starting from the valence wave function. The LF- 
conjugated operator, II(K) is given by 

7T(K) :=g-\K) \G L {K) (26) 

with G L (K) := G a (K) + G (K)V(K)g(K) = G (K)G- 1 (K)G{K) , that allows the "|" 
operation on the left. The reverse LF projectors make compact the relations between 
operators and states living onto a 3D hypersurface and the full 4D counterparts. For 
instance the relation between the BS amplitude and the valence component can be 
written as I*') = IJ{K)\$) and ($>\ = (<P| n(K). These relations can be applied to both 
two-boson |13II14| and two-fermion systems |20H21j with the proper choice of G, Go an d 
the on-minus-shell Go(K). Reversely, the valence component of the LF wave function 
can be obtained directly from the BS amplitude by using Eqs. H24[) and (|25|l 13], 14 ,20 

ED 

\G {K)G \K) |#) = \Go{K)G \K)G{K)G^\K)G Q {K)\ g-\K) |£) = 

= I*) • (27) 

To conclude this section, it is worth noting that the 3D valence LF wave functions are 
solutions of the squared mass eigenvalue equation: 

g-\K) |*> = , (28) 

with suitable boundary conditions for bound and scattering states, respectively. In 
particular, the LF scattering state is the solution of the inhomogeneous equation, 

K) = \t>o)+9o(K)w(K)\<P+), (29) 

with outgoing boundary condition. 
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4 Hierarchy equations for LF Green's functions 

In order to gain deep insights in the Fock content of the dynamical quantities involved 
in the approach presented in the previous Sections, one can combine the investigation 
of the Fock structure of a Hamiltonian (actually a hadronic Hamiltonian) performed in 
Refs. |10llll| . with the LF projection technique [15] • For that purpose, the free resolvent 
is rewritten as 



ga(K) = \G {K)\:= / dk'fdk, (fc'f | G Q (K) \k^) = id 



(30) 



where the phase space operator is conveniently defined by i? := y k^(K+ — k^). It 

should be pointed out that fl makes the LF projection " |" invariant with respect to the 
kinematical subgroup of the Poincare group. For spinless particles, the free two-body 
LF Green's function is a particular case of the N-body LF Green's function given by 



9 { N) (K) 



N 



l[e(k+)9(K + 

3=1 



K 



(N)- 



+ ie 



(31) 



where K, 



(AT)- 



is the free LF Hamiltonian. 



^/j — 1 1X1 jon 

Let us consider a two-boson system. By introducing the operator O, the interacting 
LF Green's function, Eq. (|17[) . can be rewritten as follows 



( K ) = (A') + g^ (K)u(K)g^ (K) 



(32) 



where g (2) = -iQg{K)Q, v(K) = if2 1 w(K)Q 1 . From Eq. US]), the leading and 
next-to-leading order contributions to v(K) are given by 



n go (K) 



\G (K)V(K)Go(K)\ g (K)n 



(33) 



n go (K) 



\Go(K)V(K)Go(K)V(K)G (K)\ 



f2g (K) \Gq{K)V{K)Gq{K)V{K)Gq{K)\ g {K)Q 



g Q (K)a 



(34) 



Notice that u^\K) is two-body irreducible, due to the subtraction of the last term 
in Eq. (|34p . The content of the operator v^ n ' in the LF Fock-space can be investi- 
gated, within a ladder approximation, in a Yukawa bosonic Lagrangian model, £f — 

9s4>\4>l a + <7s02^ ,2<T w ith 4>ly 4>2 an d v bosonic fields. Then, the interaction vertex 
operator, acting between Fock states differing by one quantum a, has matrix element 
given by e.g., 



(qk a \v\k) = -2(2tt) 3 5 3 (5 + k a - k) 



9S 



Vg+k+k+ 



(35) 



where the LF momenta are indicated by the convention: q = {q + ,ci±}- The states 
are normalized according to (0. The effective interaction v{K), up to next-to-leading 
order in v can be obtained from the suitable LF-time ordered diagrams, since the LF 



10 



projections, "|", allows one to play a game analogous to the case of the non relativistic 
perturbation theory in the Fock space. Then, one gets [lOUlllfTS] 

v(K) » v<- 2 \K) + v^(K) = vgf\K)v + vg {3) {K)vg^ (K)vg {3) {K)v . (36) 

By reminding that in the Yukawa model, the BSE kernel V(K) contains two interac- 
tion vertexes, v, according to Eqs. (|33[) and Eq. (|34[) one can perform the following 
identifications: i) v {2) = vg^ 3) (K)v and ii) i/ 4) = vg^ 3) {K)vg^ ] ' (K)vg {3) {K)v. Such 
expressions straightforwardly show the Fock content of each term. In particular, the 
presence of <?q 3 ' and g^ points to an intermediate propagation of three and four 
bosons, respectively. Once the previous analysis is performed at any order, it turns out 
that in general 

v(K) = vg {3] (K)v . 

In turn, from Eq. (|36|l . one sees that 5 ( 3 ) is coupled to the four-body Green's function, 
which should be coupled to the five-body one, and so on. By an obvious generalization, 
one can construct a hierarchy of coupled equations for the interacting bosonic Green's 
functions, viz 

gV\K) = gW(K)+gW(K)vgW(K)vg( 2 \K) . . . 

g w {K) = g w {K) + ff w { kW n +v {K)vg m {K) (37) 

Those coupled equations encode the full Fock-space content of the QP expansion, 
within the LF projection framework. 

An analogous study can be carried out for the two-fermion system, by adopting 
the Yukawa model given by £.f = gs^'i^icr + gs^l^l a - with ty\ and $2 being the 
fermionic fields. The interaction vertex operator, acting between Fock states differing 
by zero, one and two it's, has matrix elements given by 

{(q,s')k<r\v\(k,s)) = -2m(2iv) 3 S 3 (q + k <T -k)—. 9S 0{k+)u(q, s')u(k, s) (38) 

^q+k+k+ 

((q, s')k'M{k, s)k a ) = -2(2^) 3 <5 3 (g + k' a -k- k a )S s , s e(fc f )fl ffi (39) 

((q, s')k' a k„\v\(k, a)) = -2{2-*) 3 5 3 {q + k' a + k a ~ k)5 s , 3 — £g= g(jg)gj*£) (4Q) 

yjk'+k+ k+-k+ 

for fermion states properly normalized. The instantaneous terms in the two-fermion 
propagator give origin to Eqs. (|3U1) and (gUJ). Since v{K) has terms that couple sectors 
of the Fock space that differ at most by two sigma's [14], then one gets the following 
expression for the coupled set of Green's functions 

g {2) (K) = g {2) (K) + g {2) (K)v [</ 3 > (K) + g^ (K) + g& (K)vg^ (K) 
+ g {i \K)vg^(K)\ [ vg^(K) , ... 

9 (N \K) = gW(K)+ g W(K)v (K) + fl ^+ 2 ) (K) + g <™> (K)vg^ (K) 



+ g^ N+2 \K)vg^ N+1 \K)\ vg^(K) , . . . (41) 
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It is important to notice that truncating in the Fock space the effective interaction 
v is different from truncating the coupled set of Eqs. (|37|l or (|41[) . This can be easily un- 
derstood by considering the two-boson case and restricting the intermediate-state prop- 
agation up to four-particles, namely retaining up to g^ (K) . Then, one gets g {2) (K) ~ 
fff {K)+g ( ^ (K)vgW (K)vgW {K ) with g^ (K) ~ 5 <j 3) (K)+g® (K)vg^ (K)vg^ (K), 
where one has propagations up to four particles, given the presence of g {3) (K). On the 
other side, from Eq. (|36p . one has vg^ (K)v, without higher Fock propagations. 



5 LF projection of three-body BS equations 

The approach briefly revised in Sects. [2] and [3] has been extended to three-particle 
systems in Ref. [19]. In the last years, within a LF framework, 3-body systems have 
been i) investigated within zero-range models |35II36| . ii) applied to the description of 
the nucleon [37] and iii) adopted for analyzing the quark mass effects in heavy baryons 
[38| . It is very important to notice that the QP expansion allows one to systematically 
deal with higher Fock-state contributions to the dynamics of the three-body valence 
component, in fully analogy with the two-body case. 

The starting point of the investigation performed in Ref. [T5] is the three-body 
Bethe-Salpeter equation for the transition-matrix, within the ladder approximation, 
viz 

T = V + VG Q T ■ V = ^Vi ; Vi = V$Sf l , (42) 

(2) 

where V , stands for an interaction between particles j and k corresponding to 2-body 
irreducible diagrams. Si is the individual particle propagator and Qq = S1S2S3. The 
complexities produced by the spin degrees of freedom are omitted in what follows. 

The same QP formalism shown in Eqs. (|11|) and (|12f) can be applied to the three- 
body BSE, but with the obvious extension to the three-body case of the free Green's 
function, namely Qq = S\S% — > Go — S1S2S3. In order to get the three-body LF 
Green's function go(K), one has to project onto the LF hyperplane by integrating 
over two independent minus momentum components, i.e. Go := Go||<7q~ 1 ||Go with 
go '■= 1 1 Go 1 1 and Go the generalization to the three-particle case of the on-minus-shell 
two-body propagators (see Eq. (0 for bosons and Eq. (O for fermions). The double 
bar "|| "operation on the right or on the left means 

" G ° := J dk i dk 2 ( k i ' k 2 \ G °' G °H := f d Pl d P2 G br>i>2 )■ ( 43 ) 

The Faddeev decomposition of the QP as W\ = Vi + V{Aq y\ Wj, leads to the 

components of the effective interaction Wi := rioWillo := 3o~ 1 ||GoW / iGo||3o~ 1 - The 
free reverse LF time operators for the three-body system, 77o and 77 0, are introduced 
in analogy to the two-body case, with the difference that they now contain a double 
integration. Finally three-body LF transition matrix reads as follows 

3 3 

t := 77 T77o = ^ m » + $Z Wl9ot ' (44) 

i=l i=l 

where the standard Faddeev decomposition, t = leads to a couple set of 

equations t, = t< + hg^tj + t k ) with t» = (1 — w l g )~ 1 w i . 
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One is tempted to identify tj with the two-body subsystem 3D transition matrix, 
however this is not possible. It contains irreducible three-body terms from the point 
of view of the LF global propagation. To clarify this point, let us consider an example 
with the interaction Lagrangian density: £j = ^- =1 g4>\4>iCJ, where three different 
spin zero bosons exchange a scalar quantum a Moreover, let us simplify our discussion 
assuming that Vi corresponds to the 4D ladder one-boson exchange. The expansion of 
Wi, e.g., up to next-to-leading order is given by 

Wl = n Vjlo + n ViA Vjlo + IIoViAoiVj + V k )7T + ■■■ , (45) 

The leading-order 77ol/;-/7o corresponds to a LF two-particle interaction in the presence 
of the i-th boson, and it is built through the coupling of three- and four-particle Fock 
sectors. Pictorially, one has an intermediate propagation of a a four-particle state, 
i.e. three 4>i bosons and an exchanged quantum a, between an initial and final three- 
boson free propagations. The term noViAoVillo corresponds to two-body stretched 
boxes involving two boson out three, with the third one, the i-th boson, acts as a 
spectator. The term IJoViAoVjllo corresponds to an induced three-body force due to 
the elimination of the relative LF-time between the particles. It should be pointed out 
that the induced three-body forces, from the 4D point of view, are quite different from 
the intrinsic three-body forces, and play a very important role in the determination 
of the three-body dynamics, as shown by Karmanov and Maris in [39], where the 
calculation of the three-boson bound-state mass has been presented. 

6 Electromagnetic Current and LF Ward-Takahashi Identity 

The electromagnetic current operator plays a central role for the phenomenology, and 
therefore it deserves a detailed analysis within any relativistic framework. In particular, 
one should pay attention to the fulfillment of the Ward-Takahashi identity (WTI) [3D] , 
as well. For a system of two charged particles 1 and 2, in the Minkowski space, the 
WTI reads as follows 

Q»J"(Q)= [£T\ei] + (l<->2), (46) 

where the charge operator for particle i has matrix elements given by (ki[ei\pi) — 
eiS 4 ^ (ki — Pi — Q) ■ The full Green's function of the interacting two-particle system is 
a solution of G(K) = Qo(K) + Qq(K)V (K)G(K) ■ The current operator may contain 
two terms, a free contribution and an interacting part, viz 

J"(Q) = Jf{Q) + J?{Q), (47) 

It is worth noting that the free term, Jq(Q), leads to the impulse approximation 
derived by Mandelstam [41], where self-energy insertions/vertex corrections were dis- 
regarded. 

Once the relation between the BS amplitude and the 3D LF valence component (see 
Sect. [3]l has been established, the LF em current operator can be constructed from the 
matrix element of the 4D current (see |20II21| ). For both scattering and bound states 
one has 

J 11 iff) Wi) = (0/1^(^,^)1^) , (48) 
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where ^{K f , Ki) is the 3D LF current, that acts on the valence wave functions and 
includes two-body operators. It is given by 

f{K f ,Ki) ^HiKWWlIiKi) . (49) 

A 3D LF electromagnetic current operator for two-boson interacting systems, acting 
on the valence state and fulfilling a WTI, has been also obtained by using the method 
of gauging equations [42]. The achieved result is in agreement with Eq. (1491) . 

In order to implement the WTI for the 3D current one starts with the corresponding 
relation for the 4D current, Eq. (|46|l . and then one applies the reverse LF projector, as 
suggested by Eq. (J49j). The WTI satisfied by j^{Kf,Ki) is given by 

Q tl f(K f ,K l )=T7(K f ) [G-^K^-eG-^Ki)] n(K t ) , (50) 

with tt= e\ + e2. After applying some formal manipulations, illustrated in great detail 
in [2T], one gets 

Q^(K f ,K t ) = g-\K f ) Q L LF - Q R LF g' 1 ^) (51) 

where there are the fully interacting 3D LF Green's functions, labeled by the total initial 
and final momenta, and the left and right LF charge operators have been introduced, 
according to the following definition 

Qlf = \G (K f ) Go\K f ) en(Ki) = \G (K f ) G^ 1 (K f ) e n (K z ) , (52) 

Qlf =n(K f )eg^ 1 (K l ) GoiK^^TfoiK^eg^iK,) G (K t )\ . (53) 

with the same formal assumption indicated below Eq. ([8|. In the case of fermions, the 
explicit expression, e.g. for particle 1, is given by 

QlLF = A +( k lon)^Tft ^XLFA+(ki on )A + {k2on) , (54) 
k l 

QlLF = A+(ki on ) e 1LF ^-3-7 1 f A + (ki on )A + (k 2o n) , (55) 

where the notation eiLp indicates the 3D LF counterpart of the 4D operator e, with 
matrix elements given by 

(fc;+,k' 1± |e liF |fc+k lj± ) := ei 5 3 (fci -k x -Q) . (56) 

The normalization condition can be obtained sandwiching the operator j + m/k + be- 
tween LF spinors. 

In the case of a spin-zero boson system, one simply has Qxlf = Qxlf ~ ^ILF- 

It is very important to stress that the LF charge operator are interaction free, 
within the approach of Refs. |20H21j . where no self-energy corrections were included. 

From Eq. (|5ip , current conservation straightforwardly follows by taking the matrix 
elements between 3D LF valence components, solutions of the wave equation (|28[) and 
noting that the left and right charge operators do not contain any le dependence. 

By multiplying both the left and right hand sides of Eq. (|51[) by g(Kj) and g(Ki), 
respectively, one gets 



Q»g{K f )f{K } ,Ki)g{Ki) = QLF9(Ki)-9(K f )Q% F , 



(57) 
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which corresponds to the LF projection of the 5-point function without instantaneous 
terms in the external legs in the case of fermions. For bosons this identity was also 
derived in |42| . 

Within a field theoretical approach, the Poincare covariance of is ensured. There- 
fore, the covariance properties of the matrix elements of the operator j^, are fulfilled, 
since all the matrix elements of the lhs of Eq. (|48p are properly related through the 
Lorentz transformations. But for a truncated QP expansion, see Eq. {TS}, it is expected 
that the full covariance of the description will be lost, at some extent. A quantitative 
analysis of such a Poincare covariance breaking can be systematically carried out by 
considering more and more terms in the QP expansion. In any case, even after introduc- 
ing a proper truncation, the corresponding LF current fulfills a WTI, where truncated 
LF Green's functions appear [20U21] , The corresponding current conservation can be 
retrieved by using the valence wave functions obtained from eigenequations with the 
inverse of the truncated LF Green's function. Although the explicit expression of the 
truncated LF current is rather involved, the physical picture that arises is quite sen- 
sible: the truncated LF two-body current, that fulfills the proper WTI, is generated 
by attaching the photon in all the possible ways to the truncated effective interaction 
operators, that are present in the truncated reverse LF projectors (see Eq. (dU) and 
are irreducible with respect to LF two-body propagation. 

Within a framework kinematically Poincare invariant, the rich phenomenology of 
the electromagnetic processes can be addressed by adopting the approach presented 
in this Section. It is interesting just to remind few issues. In the fermionic case, the 
instantaneous contributions play an essential role at the formal level, and therefore 
one should be eager to investigate possible signatures of those terms. Furthermore, for 
both bosonic and fermionic systems, the choice of a frame different from the celebrated 
Drell-Yan one, namely a frame where the plus component of the momentum transfer 
is not vanishing, allows one to study the Fock content of the system state, by coupling 
small components to large components (for an analysis of the physical impact on actual 
cases in Hadronic Physics see, e.g., Refs. [17| . [43] and [44] ) . 

7 Summary and perspectives 

In this review, we have given a short presentation of the LF projection method, based 
on a combination of the integration on the minus component of the constituent four- 
momenta and the quasi-potential formalism [12] ■ This approach allows one to formally 
establish a one-to-one relation between the Bethe-Salpeter amplitude and the 3D LF 
valence component of a bosonic or fermionic system, and it makes natural to address 
the issue of the Fock content of the dynamics governing the system under investigation. 
In particular, in Sect. [4] it has been shown that by using the quasi-potential approach, 
in the spirit of the "iterated resolvents" suggested by H.C. Pauli [101111] , one can 
construct a set of coupled equations for the LF resolvents, that allows one to arrange 
a sort of tomography in the Fock space of the Bethe-Salpeter equation. Finally, the 
extension of the LF projection method to the three-particle case has been outlined. 

From the phenomenological point of view, the last Section contained the most 
relevant topic. There, we briefly presented the derivation of the conserved 3D LF elec- 
tromagnetic current operator and the associated Ward-Takahashi identity. The issue 
of the truncation in the Fock space of the LF current, and the consequent elaboration 
for obtaining a truncated Ward-Takahashi identity has been discussed. In particular, 



15 



it has been pointed out that in order to fulfill the current conservation, the valence 
wave functions have to be solutions of the suitable eigenequation, obtained from the 
truncated Green's function. 

In conclusion, we would mention that, within the framework of the LF projection 
method, the calculation of the deuteron electromagnetic form factors, with two-body 
currents obtained from a pion exchange interaction, and the investigation of the re- 
lation between the analytic properties of the Nakanishi representation and the Fock 
decomposition of the BS amplitude (see e.g. [27II28| ) are in progress. 
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